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Abstract 


This pilot study was to assess the high temperature and zero baryon density region of quantum 


chromodynamics (QCD) phase diagram with thermal photon emission, where the nature of QCD phase transition is 


ambiguous. Based on a (3+1)-D ideal hydrodynamical model to describe macroscopically the collision system, 


thermal photons emitted from Pb+Pb collisions at 2.76 TeV are investigated. The result reveals that photons from 


heavy ion collisions at high energy and centrality are possible to distinguish the structure of the hot dense matter, in 


QGP phase or hadronic phase, thus may provide an approach to explore the nature of this finite-temperature QCD 


transition (that is, first-order, second-order or analytic crossover). 
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1 Introduction 


Quantum chromodynamics (QCD) is the theory of the 
strong interaction, explaining (for example) the 
binding of three almost massless quarks into a much 
heavier proton or neutron. The standard model of 
particle physics predicts a QCD-related transition. At 
low temperatures, the dominant degrees of freedom 
are colorless bound states of hadrons (such as protons 
and pions). Because of asymptotically free, at high 
energies or temperatures, hadrons break up into quark 
degrees of freedom. Despite enormous theoretical 
efforts, the nature of this finite-temperature QCD 
transition (that is, first-order, second-order or analytic 
crossover) remains ambiguous. Here we explore the 
nature of the QCD transition with thermal photon 
emission in heavy ion collisions. 

Photon emission rates from different phases, 
hadronic gas and quark gluon plasma, are investigated. 
According to Kapusta’s pioneer investigation''!, the 
conclusion is that photon emission rates from the two 
different phases are comparable at temperature around 
200 MeV. However, in QCD phase diagram, the 
temperature involving the finite-temperature QCD 
transition has a quite wide range, from T, ~170 MeV 


to a few times of T.. Therefore photon emission rates 
from two phases at temperate ranged from 200 MeV to 
much higher temperature, i.e., 700 MeV, have been 
investigated. On one hand, it is hard to accept the 
existence of hadron gas at temperate higher than Te. 
On the other hand, the matter situation is not a pure 
QGP phase according to the most popular lattice 
results, a cross-over phase transition. Therefore, we 
still consider the part other than QGP matter as a 
hadronic matter and take photon emission rate from 
hadronic gas to estimate the unknown emission rate. 
The temperature-dependent photon emission 
rate study tells that at very high temperature, the 
photon emission rates from different phases differ 
from each other with a factor of several magnitudes. 
This will make a clear identification of the nature of 
this finite-temperature QCD transition possible. 
However, one needs a real measurement based 
on a certain experiment to distinguish the different 
options of this finite-temperature QCD transition. 
Therefore, a high temperature matter should be created, 
and the photon emission should be measured. The best 
choice is Pb+Pb collisions at 2.76 TeV, which makes 
the existence of high temperature matter possible. 
Then based on a (3+1)-D ideal hydrodynamical model 
to describe macroscopically the collision system, 
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thermal photons emitted from Pb+Pb collisions at 
LHC energy 2.76 TeV with two options phase 
transitions are investigated. 


2 Approach and results 


In this section we will introduce our approach in detail. 


Here we need several segments of this study. We 
introduce them one by one in the following. 


2.1 Thermal photon emission rates 


Thermal photon production is obtained by integrating 
the photon emission rate R (number of reactions per 
unit time per unit volume which produce a photon) 
over the space-time history of the expanding hot and 
dense matter. In this section we study the photon 
emission rates from different phases of the hot dense 
matter. 

The spectral photon emissivity directly reflects 
the dynamics of real photon production reactions in 
thermalized matter. Commonly employed formalisms 
are finite-temperature field theory and kinetic theory. 
As systematically studied by Kapusta et al., the 
thermal emission rate of photons with energy E and 
momentum p from a small system (compared to the 
photon’s mean free path) is 


dR -2 
dp (2x) 


Ryu 1 


“ exp(E/T)-1? 


(1) 


where TA is the retarded photon self-energy at 
finite temperature T. This formula has been derived 
both perturbatively and nonperturbatively. It is valid to 
all orders in the strong interaction. If the photon self- 
energy is approximated by carrying out a loop 
expansion to some finite order, then the formulation of 
Eq.(1) is equivalent to relativistic kinetic theory, 
where the emission rate of photons with energy E and 


momentum p from a process of type 1+2—3+y 


reads as 


E p, 
ajm e e 
d p iq (27y 2E, 
x5 (p} + pt — py — p”) 


2 fA) ALAM + f(E) 
2(22)° , 


(2) 


x|M| 


where f’s are the Fermi-Dirac or Bose-Einstein 
Eq.(2) is 
convenient if the scattering amplitude, M, is evaluated 


distribution functions as appropriate. 
in a perturbative expansion. Non-perturbative (model) 
calculations at low and intermediate energies, on the 
other hand, are more amenable to the correlator 
formulation, Eq.(1). In the hadronic medium, e.g., 
Il, can be directly related to vector meson spectral 
functions within the vector dominance model (VDM). 
Instructive investigation on photon emission rates 
from both QGP and HG phases can be found in 
Refs.[1,2]. 

The thermal rate from a quark-gluon plasma is 
computed using the kinetic theory formalism for the 
simplest two-to-two scattering diagrams such as the 
QCD Compton process gg—yq and annihilation 
qq > gy”. As noted in Ref.[3], this does not yet 
comprise the full result to leading order in the strong 
coupling constant @,. Due to collinear singularities, 
bremsstrahlung as well as pair annihilation graphs 
contribute at the same order as the resumed 2—2 
processes. The full result, which also necessitates the 
incorporation of Landau-Pomeranchuk-Migdal (LPM) 
interference effects, has been computed in Ref.[4,5] as 


n 2 
pe _< G Ad, q 1 
d p mle] 2x e +1 


(3) 


fo) Hn(2x) H Ca (xX) + Crrems (X) + Conn (20) 


of the three 
functions C. The corresponding results at different 


with convenient parameterizations 


temperatures are illustrated in Fig.1 as solid lines. 
Photons can also be produced in a hadronic 
phase, from several elementary interactions. The 


12 n 
[L2] comes from the reactions 


dominant contribution 
mm—py and mp—-ay. The decay p—n+n+y also 
contributes significantly. Interactions involving 
strange mesons or baryons can also produce photons, 
but these contributions are relatively small because of 


the phase-space suppression due to their big masses. 
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Fig.1 (Color online) Photon emission rates at temperatures 7=200 MeV, 300 MeV, +*+, 700 MeV from QGP phase (solid lines) and 


hadronic phase (dotted lines). 


The situation of thermal photon radiation rates 
from a hadronic gas is uncertain, due to difficulties 
related to the strong coupling and the masses of 
hadrons. The study is usually carried out within 
effective Lagrangians. Constraints on the interaction 
vertices can, to a certain extent, be imposed by 
symmetry principles (e.g, e.m. gauge and chiral 
invariance). Coupling constants are estimated by 


adjusting to measured decay branchings in the vacuum. 


Thus, for the temperature ranges relevant to practical 


applications, the predicted emission rates are 
inevitably beset with significant uncertainties, and 
therefore a careful judgment of the latter becomes 
mandatory. We will use the results of the MYM 
calculation! and plot the corresponding parameterized 
rates in Fig.1 as dotted lines. 

We can see that at temperature below 
500 MeV, the emission rates from the two different 
phases are comparable. But the temperature reaches 
600 MeV or even higher, the emission rate from 
hadronic gas is several magnitudes higher than the rate 
from QGP phase. This 


identification of what is the matter at the relevant 


makes possible an 
region of the QCD phase diagram, thus to tell the 
order of phase transition. 


2.2 Evolution of a heavy ion collision system 


The expanding local-thermalized matter created in 
heavy-ion collisions, is treated by employing three- 


7) via the flow velocity u, 


dimensional hydrodynamics 
the energy density £, pressure P, the entropy density s 


and the baryon number density n, as functions of 


the space-time position (y, t, r, o), with 7 being the 
space-time rapidity, and r, g being the transverse 
coordinates. Those quantities are governed with a 
hydrodynamical equation: 


gT” =0, (4) 


where energy-momentum tensor can be decomposed 
as T” =(€+ P)u“u" — Pg” in ideal hydrodynamics. 
This hydrodynamical equation describes the collective 
motion of the collision system from an initial time 
until a freeze-out condition. 

An event generator EPOS has been used to 
Additionally, a 
dynamical equation, i.e., the relation between energy 


construct the initial condition. 
density and the pressure is needed to close the 
hydrodynamical equation. The equation of state used 
is Ref.[8]. For more details of the solution of hydro 
equation, one can read Ref.[7]. In Fig.2, the resulted 
energy density at the initial time t=0.35 fm/c from 
Pb+Pb collisions at Jsw =2.76 TeV at centrality 
0%-—40% (Fig.2a) and 40%-80% (Fig.2b) are shown. 
In this model, centrality is defined with impact 
parameter b, and constrained with experimental data 
such as rapidity distribution of produced hadrons and 
multiplicity distribution function. It is equivalent to the 
definition based on Glauber model. Here the impact 
parameter is be[0,9.85) and [9.85,13.93) fm for 
0%-40% and 40%-80%, respectively. One can see 
more energy has been deposited in the midrapidity 
region to form the hot dense matter in more central 


collisions. 


050524-3 


202306.00568v1 


chinaXiv 


ChinaXiva ERAT! 


LIU Fuming / Nuclear Science and Techniques 24 (2013) 050524 


(a) 


y / fm 
(= 
e / (GeV/fm°) 


aS « = oa A 
x | fm 


8 
(b) 
6 35 
4 30 
2 2am 
E 
- 
Ey 2 = 
> J 
2 15 Pa 
-4 10 
-6 5 


= @g a 2° bt a A Ue: CS 
x / fm 


Fig.2 (Color online) Obtained energy density for Pb+Pb collisions at [s 7276 TeV at the initial time t)=0.35 fm/c, centrality 


0% 40% (a) and 40% -80% (b). 


To calculate particle production, another 
equation of state, i.e., the relation between energy 
density and temperature is needed. For hadron 
production, this relation is only used at the Freeze-out 
condition. But for thermal photon emission, this 
relation is used through the whole hydrodynamical 
evolution. The relation between energy density and 
temperature is obtained based on the structure of the 
hot dense matter, i.e., a QGP or a hadronic gas, or 


anything. In this work, we still take the relation 
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between energy density and temperature from lattice 
calculation”. The obtained temperature for Pb+Pb 
collisions at „fs, =2.76 TeV at the initial time 
T=0.35 fm/c is shown at Fig.3 for centrality 0% 40% 
and 40%-80% respectively. We can see temperature 
as higher as 900 MeV has been reached in 0%—40% 
centrality, however, only in a small volume and for a 
short time interval. At most space-time points, the 
temperature will be lower than this peak value. 
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Fig.3 (Color online) Obtained temperature for Pb+Pb collisions at [sn 2-76 TeV at the initial time t9=0.35 fm/c, centrality 0%- 


40% (a) and 40%-80% (b). 


2.3 Thermal photon production 


Thermal photon production is obtained by integrating 
the photon emission rate R (number of reactions per 
unit time per unit volume which produce a photon) 
over the space-time history of the expanding hot and 


dense matter: 


dN 
dyd’ p, 


is the four-momentum of produced 


= | d*xR(p"u,,T), (5) 


where p“ 
photon. The flow velocity u“ and the temperature T 
are based on the above calculation. The photon 
emission rate covers emission from both QGP phase 
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and hadronic gas, 
R= foa RON’ += foc) R”, (6) 


where focp = Joop (¥) is the fraction of QGP phase 
at a given space-time point. 

There are two options of fp. Option one 
(optl) is consistent with the dynamical equation of 
state, the equation of state from lattice QCD", 
where pressure is calculated”! as p = P, +A(Py - Py): 
one obtained the temperature-dependent fop =1-/. 
This makes hadronic gas appear at very high 
temperature, with very small fraction, ie, about 4% at 
700 MeV. 

Option defined 
according to first order phase transition, which was 


two (opt2) of foar is 


widely used in Ref.[9], where no fraction of hadronic 
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gas appears at temperature higher than 200 MeV. 

Due to the two options of focp, the 
transverse momentum spectra of thermal photons at 
midrapidity from Pb+Pb collisions at 2.76 TeV are 
plotted in Fig.4 for centrality 0%—40% (Fig.4a) and 
40%-80% (Fig.4b). Dotted lines are results based on 
option one of foar. Solid lines are based on option two. 
At a later time, the system temperature gets cold and 
the emission of the two options is expected 
comparable due to the emission rates from the two 
different phases. Thus photons emitted at time later 
than 7=0.35 fm/c are counted as a comparison, i.e., 
from 0.75 fm/c, 1.15 fm/c and 1.55 fm/c, plotted as 


different dots and colors of lines. 
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Fig.4 (Color online) Transverse momentum spectra of thermal photons at midrapidity from Pb+Pb collisions at 2.76 TeV. 


3 Conclusion 


The comparison of the transverse spectrum of thermal 
photons from Pb+Pb collisions at 2.76 TeV reveals 
that high collision energy and centrality are needed to 
And the early 
evolution is the most sensitive stage to distinguish the 


create high temperature matter. 


structure of the hot dense matter, in QGP phase or 
hadronic phase, thus may provide an approach to 
explore the nature of this finite-temperature QCD 
transition (that is, first-order, second-order or analytic 
crossover). 
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